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Abstract 

We present predictions for the non-singlet Drell-Yan longitudinal spin cross 
sections and spin asymmetry, All, in proton-proton collisions at large pt at 
the RHIC energy of 200 GeV at next-to-leading order QCD. The higher order 
corrections to the non-singlet polarized cross section, are sizeable and 

similar to those found for the unpolarized cross section. The non-singlet asym- 
metry parameter, A^f, is very stable against higher order corrections and is a 
direct measurement of the non-singlet (i.e., valence) polarized quark distribu- 
tions within the proton. 
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1 Introduction 



The study of the spin content of the nucleon at hadron/hadron colliders is one of the 
most interesting aspects of QCD which is increasing receiving attention both from the 
theoretical and from the experimental side. 

Among the most interesting questions that a p-p collider can help us answer is the 
issue of the distribution of the proton spin among its constituents and the measure- 
ment of various polarized parton distributions. In this new setting, the presence of 
a direct gluon coupling, absent in polarized Deep Inelastic scattering processes, will 
allow a direct test of the gluon contribution to the polarization of the nucleon, with 
measurements covering a wide x (Bjorken-x) range. We mention that, especially at 
smaller x values, the gluon polarization AG is strongly model-dependent and, there- 
fore, poorly known from DIS experiments alone. 

Along these lines, in the last few years, several studies regarding the analysis of 
the scaling violations in polarized p-p collisions have been presented, and considerable 
progress has been made in the theoretical investigation of various important polar- 
ized processes. For instance, it has become now clear that inclusive prompt photon 
production [I]] (p p — > 7 + X), photon-plus-jet production (jp p — > 7 + J) ||, and 
polarized production of 2 jets (p p — > J J) are among the most important processes 
to be studied at RHIC, the Relativistic Heavy Ion Collider at Brookhaven. 

The corresponding cross sections show sizable asymmetries (see for instance 0) 
with good hope for the experimental program. We mention also that, beside the usual 
asymmetries for total cross sections or large px cross sections, other observables of 
the final state can be used as a way to test the size and the sign of AG. 

In this respect, the study of the event-structure of the final state and the asymme- 
tries for various correlations (for instance rapidity correlations, angular correlations in 
the photon/jet system) have been shown also to be a source of important information 
on the polarized parton distributions and can help us discriminate among the various 
different models M. 

More recently, lepton pair production via the Drell-Yan mechanism has also re- 
ceived considerable attention. The process has been historically a clean way to test 
QCD and the parton model and continues to be so also in the case of polarized col- 
lisions. In a seminal paper HJ], Ralston and Soper pointed out that the factorization 
formula for parton distributions for Drell-Yan includes, beside the usual longitudi- 
nally polarized parton distributions, also parton distributions for transverse spin. In 
particular the study of the transverse spin distribution is one of the challenging tasks 
that the STAR and PHENIX collaborations at RHIC are going to undertake within 
the next few years. 

On the technical side, we remind that although formal proofs of factorization are 
still missing, direct calculations show that the Drell-Yan cross sections does factorize, 
as shown in nontrivial cases (for instance for the nonzero pr distributions ||). This 
opens the way to an accurate study of this important cross section both for longitu- 
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dinal and transverse spin. At pr = 0, where px describes the transverse momentum 
of the lepton pair, the calculations simplify drastically and various results have been 
reported in the literature (see for instance [7L M for the longitudinal case). 

In previous works || , we have analyzed the spin dependence of hard scattering in 
Drell-Yan (with longitudinal polarization) at parton level to NLO QCD, at nonzero pt 
of the lepton pair. In this work we proceed with the analysis at hadron level of these 
former results, thereby providing predictions for the behavior of the cross section at 
RHIC energies. Since our NLO analysis is limited to the non singlet sector, our study 
is exploratory in nature, and we have limited our attention to just one particular 
center of mass scattering energy (200 GeV), right in the middle of the planned RHIC 
energy range. We will also limit our interest to a single set of parton distributions, 
specified below. Our interest, in this work, is limited to assess the stability of the full 
0(a 2 ) calculation, predict the asymmetries of various helicity dependent amplitudes 
and show that the corresponding K factors are comparable to the unpolarized ones. 

The asymmetries turn out to be sizable and can be larger than what estimated 
here once the gluon contribution is fully included. We hope to return to a re-analysis 
of the complete process in the near future. 



2 Structure of the Differential Cross Section 

We begin by examining the contribution to the production of virtual photons with 
invariant mass, Q, in hadron-hadron collisions, 

Ht + H^^ + X, (1) 

from the 2 — > 3 parton subprocess, 

a + 6->7* + c + d. (2) 

Here Hi and H 2 are incoming hadrons with 4-momenta, P\ and P 2 , respectively, 
and q is the 4-momentum of the virtual photon, 7*, as shown in Fig. 1. The 4- 
momenta of the incoming two parton a and b are labeled by p\ and p 2 respectively, 
and the outgoing 4-momenta of partons c and d are labeled by k 2 and k%. The virtual 
photon 4-momenta is given by ki = q. Conservation of energy and momentum implies 
that 

Pi +P2 = h + k 2 + k 3 . (3) 
The hadron-hadron process (|1|) is described in terms of the invariants, 

S=(P 1 + P 2 ) 2 , T=(P 1 -q) 2 , U={P 2 -qf. (4) 

In addition, it is convenient to define the scaled variables 

x T = 2q T /W and r = Q 2 /S, (5) 
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where W = y/S is the hadron-hadron center-of-mass energy and qr is the transverse 
momentum of the virtual photon with invariant mass, Q. It is also useful to define 
the following two quantities, 

X! = = ieVx§, + 4r, x 2 = = + 4r, (6) 

where y is the rapidity of the virtual photon. The 2 — > 3 parton subprocess (^) is 
described in terms of the invariants, 

s = (Pi + P2) 2 , * = (pi - <?) 2 , u = (p 2 - q) 2 , (7) 

and 

s 2 = s 23 = (k 2 + h) 2 , (8) 
where momentum and energy conservation ([3D implies 

s + t + u = Q 2 + s 2 . (9) 

In QCD the hadronic cross section is related to the parton subprocess according 

to 

da^ + H 2 ^ ^ + X;W } q T ,y) = 

Gn.^aixi, M 2 )dxiGH 2 ^b(x 2 , M 2 )dx 2 [ (ab — > 7*cd; s,t,u)\ dtdu, (10) 

\dtdu J 

where Gn^ai^i, M 2 )dxi is the number of partons of flavor a with momentum 
fraction, x± = pi/Pi, within hadron Hi at the factorization scale M. Similarly, 
G H 2 -^b( x 2i M 2 )dx 2 is the number of partons of flavor b with momentum fraction 
x 2 = p 2 /P 2 i within hadron H 2 , at the factorization scale M. In the remainder of 
this paper we will take the factorization scale to be Q and evaluate the parton sub- 
processes at the same scale Q. Using (f|) and (jfy we see that 

s = xiX 2 S, (11) 
(t-Q 2 ) = xi(T - Q 2 ) = -xiX 2 S, (12) 
(u-Q 2 ) = x 2 {U-Q 2 ) = -x 2 xiS, (13) 

which from (|9]) implies that 

xix 2 - xix 2 - x 2 xi = r 2 - r, (14) 
where xi and x 2 are defined in (||) and 

T 2 = S 2 /S. (15) 
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It is now easy to compute the Jacobian 



dx 2 dt = — — — dxtdy, (16) 

4(xi-xi) 

which when inserted into (]T0|) and integrating over x\ and s 2 gives 

da f 1 [ s ™ ax 1 

S j 2 a ( W ^T,y) = dxi ds 2 t — r 

dq^dy Jx™ m Jo [Xi — xi) 

do 

G Hl ^ a {x x , Q 2 )G H2 ^ b (x 2 , Q 2 ) s——{ab -> ^*cd; s, t, u), (17) 

dtdu 

X X X 2 + T 2 -T 

x 2 = = , (18) 

X\ — X\ 

s™ ax = A=(t-x 1 + x 1 (1-x 2 ))S, x T n = ^—^- (19) 

1 — x 2 

The maximum value of t 2 arises when x 2 = 1 in fll4l) , while the minimum value of x\ 
occurs when r 2 = and x 2 = 1. 

For 2^2 parton subprocesses such as the Born contribution in Fig. 2 one has 

do , . _ . . do , . . 

s dtdt^ S) £ ' ^ = ^ S ~dt^ a + 7 c; s ' ^ 20 ^ 

which when inserted into (|T7D results in 

do- r 1 



where 



and 



S , °. (W,QT,y) = / . da?i 
dathdv Jx™ m 



dq\dy ' ' Jx™" (xi — x\) 

Gh^x^Q^Gh^x^Q 2 ) S ^(a6^ 7 *c; S ,t), (21) 

where in this case s 2 = and s + t + u = Q 2 . 

The Drell-Yan differential cross section for producing a muon pair with transverse 
momentum qj- and invariant mass, Q, at rapidity y in the hadron-hadron collision, 

F 1 + F 2 ^( 7 *^/iV~)+^ (22) 
at center-of-mass energy W, is given by 



S dQHq 2 dy {HlH2 ^ + X; ^ ?T ' y) 



" ^- F ^-(fr 1 fl 2 -7* + X;W, a r ) y), (23) 



\3irQ 2 J dq\dy 

where a is the electromagnetic fine structure constant and where the virtual photon 
differential cross section is given by ( |TTD or fl2"T|). 
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2.1 The Non-Singlet Cross Section 

We define the "non-singlet" cross section at the hadron level to be the antihadron- 
hadron minus the hadron-hadron cross section as follows: 

°NS = °H lH2 - VHtfa, (24) 

and at the parton level to be the quark- ant iquark minus the quark-quark cross section, 

&NS = &qq - <?qq- (25) 

Since the antihadron-hadron cross section is larger, this definition yields positive val- 
ues for <Jns- For parton distributions, on the other hand, we use the customary 
non-singlet definition, 

^H-tq 

In proton-proton collisions the non-singlet cross section is given by 



G 1 ^ — Gn^q — Gn^q- (26) 



NS _ _._ _ _ 
°pp — U VV U PP 

= (Gp^qGp^g + Gp^qG p -^q) ® &qq 
+ (Gp^qGp^g + Gp^qGp^q) <g) & q q 
— (Gp-^qGp^q + Gp^qGp-ig) ® Oqq 
(Gp^qGp—,q + G p^qG p—*q) ® &qq 

= G™ q G™ q ®a™, (27) 
where the convolution <g> is defined in ([[7]) and fl2~l"|) . 

2.2 The Spin Dependent Cross Sections 

The non-singlet proton-proton helicity dependent Drell-Yan cross section can be writ- 
ten in the form, 



G^ ] -,, M (x 1 ,Q 2 )G^ l -, q(> ,, ] (x 2 ,Q 2 ) S ^M( S , (,„,/,), (28) 
where Ai and A 2 are the helicities of the initial two hadrons, respectively, and h = 



4A1A2. In 53 the sum is over the quark flavor and the quark helicities Ai and A2, while 



the non-singlet helicity structure functions are given by 

T p[A]^q f [X] 



G p[A]-+q f [\] - Gp[A]-*q f [X] ~ Gp[A]^q f [ X }- (29) 
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For one quark flavor 



and 



(7 , 



+ (G^G^ + Gfi 4 Gfi T ) ® (31) 

where j and J, refer to helicity + | and — |, respectively, and where we have used 
simplified "convolution notation". The unpolarized non-singlet cross section is the 
average of a+l and a^ s 



a^ = K<f + <-)> (32) 



'++ iy +-> 

''ATS - 2 \ u ++ ^ U +-J 

while the non-singlet longitudinal spin difference (i.e., polarized) cross section is de- 
fined by, 



(7, 



LL 1 („NS ^NS\ /qq\ 

ns = 2 r++ ~ a +-J • ( 33 ) 



From ( p0|) and (31) it is easy to show that, for one quark flavor, 

°ns = (G™ q G™ q ) ® <r£ 5 , (34) 

and 

4s = (AG™ AG™) ® a L N % (35) 



p->q 

where the unpolarized parton distributions is the sum of the two spin states, 

and the "polarized" parton distributions is the difference of the two spin states, 

AG™ q = G™ - Gfi,. (37) 
Similarly, the parton unpolarized non-singlet cross section is the average of try f and 

°+-' 



a JV5 = 2 V a ++ + a +~) ' ( 38 ^ 

and the polarized parton non-singlet cross section is defined by, 



a NS - 2 V a ++ ~ °+-J ■ ^> 

In Ref. || we showed that to order a 2 s the Drell-Yan non-singlet spin-dependent 
parton-parton differential in (p8|) can be written as the Born term plus four (unpolar- 
ized) cross sections as follows: 

s< -J7T'{ s ^ u ^ h ) = ( X _ h ) S ( s 2) s ~lf 
dtdu at 

- ^ a lr + (1 " + (1 - - (1 + ' (40) 

arcra dtdu dtdu dtdu 
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where h = 4A1A2. Hence, to order a 2 s the Drell-Yan non-singlet unpolarized differential 
cross section for producing muon pair with transverse momentum and invariant 
mass, Q, at rapidity y in the proton-proton collision at center-of-mass energy, W, is 
given by, 

( a \ y r 1 dxi r™ ax ds2 1 

G™ q ( Xl ,Q2)G™ q (x 2 ,Q 2 ) s^(s,t,u), (41) 
S ff% t, u) = 5(s 2 )s^i + s^ + s^l + s^l- (42) 



where 



dtdu dt dtdu dtdu dtdu dtdu 

and the polarized cross section is given by, 

da LL 



,, ry , . / dx * / ds * 7 



where 



AG^faWAG^faCf) 8 ^{8,t,u), (43) 



<^o"ws / \ s daV> daf daf daV da? , . 

t— 4^ (s, t, u) = -5(s 2 )s—^ - s—^- - s— f- - s—^- - s— 44 
dtdu dt dtdu dtdu dtdu dtdu 



The unpolarized Born term, <xg and the four spin averaged cross sections, ay, ay , cf. 
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and erf, are given in Ref. [BJ and Ref.fT0[ 



3 Proton-Proton Collisions at 200 GeV 

3.1 Polarized and Unpolarized Cross Sections 

To illustrate the effect of the higher order terms in muon-pair production in hadron- 
hadron collisions, we compute the non-singlet polarized and unpolarized Drell-Yan 
cross sections in proton-proton collisions at W = 200 GeV. The non-singlet unpolar- 
ized cross section is given by (|4T|) , which in convolution notation becomes 

°%s= E (GJ^G^)®^, (45) 

q=u,d 
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where from (12) 



&ns = °b + cri + 0-2 + 03 - 0"4- (46) 
The non-singlet longitudinal (polarized) cross section given by (f43"D becomes 



LL - £ (AG^AGJ^) ® (47) 

q=u,t 



where from (HI) 



(7^5 = — (7 B - 0"! - <7 2 — 0- 3 — <7 4 . (48) 

We use GS solution A || for the polarized and unpolarized parton distributions 
(A = 200 MeV) which are shown in Fig. 3 at Q = 2 GeV. Fig. 4 shows the polarized 
and unpolarized non-singlet cross sections to leading order (Born term) and to order 
a 2 s at W = 200 GeV, Q = M = 2 GeV, and y = 0. Fig. 4 shows the "number density" 
dN / drdxxdy which is constructed from the cross section as follows: 

da ? S*-J^r, (49) 



drdxrdy 2 dQ 2 dq\dy 

using an integrated luminosity of 100/pb. 

Fig. 5 shows the non-singlet polarized and unpolarized cross sections to leading 
order and to order a 2 s at W = 200 GeV and Q = M = 2 GeV integrated over y 

da fVmax da 

dy. (50) 



drdxT Jymin drdxxdy 

Here the number density dN/drdxx corresponds to the number of events per unit 
r per unit xt in 100/pb. To compute the number of events in a bin of r of size 
0.0001 (corresponding to AM = 2 GeV at W = 200 GeV) and a bin of x T of size 0.01 
(corresponding to Apt = 1 GeV at W = 200 GeV), for example, one would multiply 
dN/drdx T by 10" 6 . 

Fig. 6 and Fig. 7 show the ratio of the full order a 2 s to the leading order cross 
sections (i.e., the "K-factors") for the polarized and unpolarized case. As has been 
previously seen, the K-factor for the unpolarized cross section is quite large. The 
polarized K-factor is almost the same as the unpolarized case. The difference comes 
only from 04 in ([46]) and ( pj|) and from the fact that the polarized and unpolarized 
parton distributions are different. However, the contribution to the cross sections 
from 04 are small, as can be seen in Fig. 8 and Fig. 9, which show the individual 
K-factors arising from the four terms in (|46| ) for the unpolarized cross section and 
from the four terms in ( f4"8"D for the polarized cross section, respectively. For both the 
polarized and unpolarized case the dominant term in the K-factor is <j\. 
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3.2 Spin Cross Sections 

The spin cross sections a+f and are calculated from a^ s and a^ L s as follows: 



n NS _ i n LL 

°++ — a NS + "iVS; 



and are shown in Fig. 10. To leading order at the parton level, 

= (leading order), (52) 
which from (^) - at the hadron level - implies an "indirect" contribution of the form 

= 2 E ( G f-T G fA + G W G fA) ® °NS (leading order), (53) 

where 

&ns = ®b (leading order). (54) 

Fig. 11 shows that the probability of finding a quark carrying helicity opposite to 
that of the proton ((jt_>j, or G[_>t) is not small especially at small x and hence at the 
hadron level a™ is not zero. 
At order a 2 s1 

a ++ — ~ 2 J2u,d (pfJi^Gf^ + <j^jG|^jJ <E> <3"4 

+2 £ u , d (Gfi T G^ + Gfi ; Gfi T ) ® (*= fl + a 4 ) , (55) 

where (jj^ is given in fl4*6|). Now there is a "direct" contribution to a+f coming from 
(T4. However, this contribution is small and the primary higher order contribution 
comes from the d\ correction to a^ s in the "indirect" contribution. 

3.3 Spin Asymmetry All 

The non-singlet spin asymmetry parameter A^l is defined by, 

NS _ NS LL 
A NS _ °++ ° + - _ a NS 
A LL — MS , ZNS ~ _£ ' 



n NS , n NS _£ ' 



and from (|4"5| ) and ( pE?D see that 



r V5 



Equation (0) implies that to leading order at the parton level, 

= and ^ns = _y (i eac ii n g or der). (5£ 
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We can make a "zero order" estimate of A?l by simply ignoring the convolution in 



and using (58) which results in 



A LL (order 0) = (A *„ Nf , ^ ; jtxttw^TWs ( 59 ) 
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)G NS G NS 4- (1)O ns G NS 



Fig. 12 shows this ratio of parton distributions at Q = M = 10 GeV versus x. The 
result is compared with the complete leading order and beyond leading order pre- 
dictions for Al£ from fl57D , integrated over y at W = 200 GeV plotted versus xt- 
At Q = 10 GeV the leading order and the beyond leading order results are almost 
identical and both resemble the zero order estimate. 

Fig. 13 shows the leading and beyond leading order results for the non-singlet 
asymmetry parameter All f° r proton-proton collisions at W — 200 GeV and Q = 
M = 2 GeV. The asymmetry is around 15% between p? of 2 and 8 GeV and there is 
very little difference between the leading and beyond leading order. As can be seen 
in Fig. 6 and Fig. 7, the K-factor for the asymmetry parameter All * s m uch smaller 
than the K-factors for unpolarized or polarized cross sections, a^ L and cr^g- 

4 Summary and Conclusions 

From (H}) and fl48|) we see that to order a 2 s at the parton level, 

-Z L ns = °ns + 26-4, (60) 

which means that except small differences due to the term 04 and small differences 
arising from the fact that the polarized parton distributions AGp_^ q in ( f4?D differ 
from the unpolarized one Gp_^ q in ([45| ), the higher order corrections to polarized cross 
section, — fr^g, and the unpolarized cross section <r^ s are identical. In both cases the 
higher order corrections are large and important. Everything that has previously been 
said about the higher order corrections to the unpolarized cross section, for example 
[10|1 , is true for the polarized case as well. 

The higher order corrections to the ratio of the polarized to the unpolarized cross 
section, 

aNS LL / £ /£-, \ 

A LL = VNslVNSi ( 61 ) 

are small since the big corrections to the numerator are canceled by big corrections 
to the denominator. As can be seen in Fig. 12 and Fig. 13, All ^ s ver y stable 
against higher order corrections. The asymmetry is around 15% between pt of 2 and 
8 GeV in proton-proton collisions at W — 200 GeV and Q = M = 2 GeV. The non- 
singlet asymmetry, All, ^ s a direct measurement of the non-singlet (i.e., valence) 
polarized quark distributions within the proton. However, to construct the non- 
singlet observables in this paper one must subtract pp spin measurements from the 
corresponding pp measurement, and it is unlikely that this will be done in the near 
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future. On the other hand, if All is measured in pp collisions at RHIC, then deviations 
from our predictions will indicate large contributions from gluons, AG p ^ g , and/or 
from the sea, AG P ^ S . This would provide the first hint as to the behavior of the 
polarized gluon parton distribution, AG p ^ g . 
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Figure 2: Radiative corrections at pt = and the lowest order contributions at 

Pt 7^0. 
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Figure 4: The non-singlet unpolarized spin-averaged (ave) cross section, a^ s , and 
longitudinally polarized (LL) cross section, cr^/g, for pp collisions at W = 200 GeV, 
Q = M = 2 GeV, and y = as a function of p?. The plot shows the "number 
density", dN/drdxTcly, in 100/pb at leading order (1) and at order a 2 s (1+2). 
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Figure 5: Same as Fig. 4 integrated over rapidity y. 
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jepjo |S|./(j9pjo puj+;s|.) 



Figure 6: Shows the ^-factors for non-singlet unpolarized spin-averaged (ave) and 
longitudinally polarized (LL) cross sections for pp collisions at W = 200 GeV, Q = 
M = 2 GeV, and y = as a function of pt in Fig. 4. The i^-factor is the ratio of the 
full order a 2 s cross section (1+2) to the leading order Born term (1). Also shown is 
the K-f actor for the non- singlet spin symmetry parameter A^l = c^vs/^ivs- 
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Figure 7: Shows the i^-factors for non-singlet unpolarized spin-averaged (ave) and 
longitudinally polarized (LL) cross sections for pp collisions at W = 200 GeV, Q = 
M = 2 GeV, integrated over rapidity y as a function of pt in Fig. 5. The i^-factor 
is the ratio of the full order a 2 s cross section (1+2) to the leading order Born term 
(1). Also shown is the K-factor for the non-singlet spin symmetry parameter A*l[ = 

rr LL Irr^ 
<J NS/ a NS- 
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Figure 8: Shows the i^-factors arising from each of the four terms in ( ^6|) for non- 
singlet unpolarized spin-averaged cross section, a^ s , for pp collisions at W = 200 GeV, 
Q — M — 2 GeV, and y = as a function of p?- The .ff -factor is the ratio of the 
order a 2 s cross section to the leading order Born term. 
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Figure 9: Shows the i^-factors arising from each of the four terms in ( f^8|) for non- 
singlet longitudinally polarized (LL) cross section, cr^-g, for pp collisions at W — 
200 GeV, Q = M = 2 GeV, and y = as a function of p?- The if-factor is the ratio 
of the order a 2 s cross section to the leading order Born term. 
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Figure 10: The non-singlet spin cross sections, a+f and cr£_, for pp collisions at 
W = 200 GeV, Q = M = 2 GeV, integrated over y as a function of pt- The plot 
shows the "number density", dN/drdxT, in 100/pb at leading order (1) and at order 
« s 2 (l+2). 
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Figure 11: Shows the ratio, G^-\/G^i, for u and d quarks at Q 
resulting from the polarized structure functions in Fig. 3. 



2 GeV versus x 
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Drell-Yan Asymmetry Parameter ALL 



0.60 x ALL =( si 9( ++ )- si 9(+-)) / (sig(++)+sig(+-)) 
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Drell-Yan Asymmetry Parameter ALL 



0.25 x ALL=(sig(++)-sig(+-))/(sig(++)+sig(+-)) 
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